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Abstract 

We consider a quantum particle in a waveguide which consists of an infinite straight Dirichlet 
strip divided by a thin semitransparent barrier on a line parallel to the walls which is modeled 
by a 5 potential. We show that if the coupling strength of the latter is modified locally, i.e., it 
reaches the same asymptotic value in both directions along the line, there is always a bound 
state below the bottom of the essential spectrum provided the effective coupling function is 
attractive in the mean. The eigenvalues and eigenfunctions, as well as the scattering matrix 
for energies above the threshold, are found numerically by the mode-matching technique. In 
particular, we discuss the rate at which the ground-state energy emerges from the continuum 
and properties of the nodal lines. Finally, we investigate a system with a modified geometry: 
an infinite cylindrical surface threaded by a homogeneous magnetic field parallel to the cylinder 
axis. The motion on the cylinder is again constrained by a semitransparent barrier imposed 
on a "seam" parallel to the axis. 



1 Introduction 

Quantum mechanics of constrained systems is experiencing a new wave of interest connected with 
the recent progress in semiconductor physics: nowadays experimentalists are able to investigate 
the behavior of electrons in structures of various shapes, at times rather elaborated. The small 
size, extreme material purity, and its crystallic structure make it possible to derive basic properties 
of these systems in a crude but useful model in which the electron is considered as a free particle 
(with an effective mass) whose motion is constrained to a prescribed subset of M d with d = 2,3, 
possibly in presence of external fields. 

On the theoretical side, this inspires questions about relations between spectral and scatter- 
ing properties of such systems and the underlying geometry and topology. A class of systems 
which attracted a particular attention are quantum waveguides, i.e. tubular regions supporting 
a Schrodinger particle. It is known that a deviation from the straight tube can induce existence 
of bound states and resonances in scattering, vortices in probability current, etc., be it bending 
[ |DE| , pEM| , pES| , |ES| , pj[|), protrusion or a similar local deformati on ]AS|, |BGRS|, |EV1| ], waveguide 
coupling by crossing ]SRW| ], or by one or several lateral windows ]ESTV| , |EV2^|EV3| (the related 
bibliography is rather extensive; the quoted papers contain many more references). 

In this paper we are going to discuss a system closely related to the last named one. It supposes 
again a double waveguide; however, the coupling between the two parallel ducts will entail now a 
tunneling through a thin semitransparent barrier rather than a window in a hard wall separating 
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them — cf. Figure |l|. To get a solvable model we describe the barrier by a 8 potential whose 
coupling strength may vary longitudinally: the Hamiltonian can be then formally written as 

H a = -A n + a(x)8(y) (1.1) 

with the barrier supported by the £-axis, where il := R x (— d 2 , di) is the double-guide strip. 

V 

di 



Figure 1: Double waveguide with a 6 barrier. 

There are several motivations to investigate a leaky-barrier waveguide pair. First of all, it is 
a generalization in a sense of earlier results, because the pierced-hard-wall case of Ref. |ESTV] 
corresponds to a = in the window and a = oo otherwise. Recall that the latter can serve 
to describe an actual quantum-wire coupler — see, e.g., |HTW , Ku| — and such a model will 



certainly become more realistic if the tunneling through the barrier of a doped semicond uctor 
material separating the two guides is taken into account. At the same time, the Hamiltonian ( |l.l| ) 
covers for various a a wide variety of situations. 

On the mathematical side, the S potential of (LI) can be treated more easily than the hard- 



wall barrier, since two operators with different functions a have the same form domain. To illustrate 
the difference, one can compare the variati onal pr oof of existence of bound states in Thm 3.1 
below with the analogous argument of Ref. | ESTV |. A deeper application of the quadratic-form 
perturbations allows us to construct the Birman-Schwinger theory for the waveguide systems in 
question, in particular, to d erive the weak-coupling behaviour of the bound states. This will be 
done in a subsequent paper |EK|. 

Let us describe briefly the contents of the paper. In the next section we shall describe the model 
and deduce its spectrum in the "unperturbed", i.e. translationally invariant case. In Section ^ 
we demonstrate that a local change of the coupling parameter will cause the existence of bound 
states provided it is negative in the mean. To ilustrate the spectral and also scattering properties 
we shall discuss then in detail the example in which the barrier function is of a "rectangular well" 
shape. In the hnal section we will show how the situation modifies if the semitransparent barrier 
is placed at the surface of a cylinder threaded by a homogeneous magnetic field. 



2 Preliminaries 
2.1 The Hamiltonian 

Let Q := R xO with O := 2 U 0\ := (— d 2 , 0) U (0, d\) be the configuration space, i.e., the part 
of R 2 occupied by the waveguide. Passing to the rational units, h — 2m — 1, we may identify 
the particle Hamiltonian H a with the Laplace op erat or away of the waveguide boundary and the 
barrier. To give meaning to the formal expression (LI) one has to specify the boundary conditions. 
At the outer edges we assume the Dirichlet condition, 

iP(x,-d 2 )=ip{x,d 1 )=0, (2.1) 

while the barrier is transversally the usual 6 potential defined conventionally as 

ip(x,0+)=ip(x,0-)=:i;{x,0), ip y (x, 0+) - ip y (x, Q-) = a(x)ip(x, 0) (2.2) 
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for any i£l — cf. | AGHH , Sec .1.3] — where the subscript denotes partial derivative with respect 
to y. The Hamiltonian domain is then 

D{H a ) := £ Wi(Q)\ V satisfies (2.1) and @ } , (2.3) 

where the function a : R — * R, assumed to be piecewise continuous, determines the shape of the 
barrier and represents the x-dependent coupling "constant" of the interaction. 

For the sake of simplicity we shall exclude the above mentioned case of a Dirichlet barrier, 
a{x) — oo at a subset of R. In that case all the operators H a have the same form domain, and the 
associated quadratic form is obtained by a simple integration by parts: 



q a [iP] := / \V^\ 2 (x,y)dxdy + / a(x)\i/>(x, 0)\ 2 dx, (2.4) 
D(q a ) := {i>€W£(Rx(-d 2 ,di))\ VxGR: ip(x,-d 2 ) = ip(x,di) = 0}. (2.5) 



The form ( |2.4| ) is obviously symmetric and it is not difficult to check that it is closed and thus indeed 
associated with H a . Hereafter we adopt the notation of [ ESTV| : d := max{di, d 2 }, D := d\ + d 2 , 
and 

minjdi, d 2 } 
v := - — . 

max{<ii, d 2 } 

Without loss of generality we may assume that d 2 < d\ = d. 
2.2 The unperturbed system 

If a(x) — a is a constant function, we can solve the Schrodinger equation H a ijj = k 2 ip by separation 
of variables. To get the transverse eigenfunctions we have to match smoothly the solutions in the 



two ducts, C 2 sui£(y + d 2 ) and C\ sin£(y — di), chosen to satisfy the condition (2J). If £di,£d 2 
are not multiples of n we get thus the following condition on eigenvalues of the transverse part of 
the Hamiltonian: 

-a = £(coted 1 +cot£d 2 ). (2.6) 



Remark 2.1 If di,d 2 are rationally related the Schrodinger equation can be also solved by £ = 
~ ^If' n e However, such wave functions are zero at y — 0, and therefore independent 

of a. In this sense they represent a trivial part of the problem. A prime example is the symmetric 
waveguide pair, d\=d 2 , where this observation concerns every solution antisymmetric w.r.t. y = 0- 
It is reasonable to concentrate on the nontrivial part only. If v = g- = |j we denote by Q w the 



subspace in L (~d 2 , d\) spanned by the solutions of Putting then H v := L (R) ® Q„ , we 



shall restrict our attention to the operator H a \ 7i v ] for the sake of simplicity we shall denote the 
restriction by the symbol H a again. The trivial part is absent, of course, if v is irrational. 

From the spectral condition ( [2.6] ) we get a sequence of eigenvalues (in the natural ascending order) 
of (the nontrivial part of) the transverse operator; we denote it as {^(oOl^Li- The corresponding 
eigenfunctions are 

Xn(y;at) = (-iyN n sm^/U^dj sin ^fu^ (y+(-iydj) (2.7) 

for y € Oj, j = 1, 2, where N n is the normalization factor chosen in such a way that x-n would be 
a unit vector in L 2 (—d 2 ,di) 1 i.e. 



K = —, . , v ~" : : : . (2.8) 

^/u^di sin y/v^d 2 + ^/u^d 2 sin ^/uildi — sin ^/u^di sin ^/u^d 2 sin ^Ju^D 



Furthermore, the Green's function of the Hamiltonian (.hi) can be written down explicitly: 



G a (x,y,x',y';k) = J2 _1_ & ik n \x-x'\ Xn{y . a) tn{y > ■ a) i (2 . 9 ) 
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where the efective momentum in the n-th transverse mode is k n := \Jk 2 — v n (a) . 

Elementary properties of the transverse eigenvalues follow from the the spectral condition ( |2.6| ) 
by means of the implicit-function theorem; we collect them in the lemma below. 

Lemma 2.2 ( a) Let {mi}^ be the sequence obtained from the set N U v^ 1 by natural ordering. 

Then ^(n— 1) < fm n _i < y/T^ < |m„ < \n holds for all n £ N\{0, 1}. 

(b) The function a i— > v n {a) is strictly increasing and continuous for all n £ N\{0}. 



3 Existence of bound states 



Depending on the choice of a, the operators (1.1) offer a variety of spectral types. In this paper 



we shall concentrate on the situation when the barrier describes a local perturbation of the system 
with separating variables considered above. The locality is at that understood as a decay of the 
function a; in other words, we shall assume that limui^oo a(x) = chq. It is important that the 
limiting value ao is the same at both directions. 

In such a case, it is easy t o local ize the essential spectrum. One employs a simple bracketing 
argument similar to that of |ESTV, Sec. II]) squeezing H a between a pair of operators with 
Dirichlet and Neumann conditions on segments perpendicular to the re-axis placed to both sides of 
the centre. By the minimax principle only the tails of the estimating operators contribute to their 
essential spectra; since the "cuts" can be chosen arbitrarily far we obtain o~ ess (H a ) — [vi{q.q), oo). 

Less trivial is the existence of discrete spectrum. It is known that any "window" in the 
impenetrable barrier induces a bound state. This fact was established first for sufficiently wide 



windows [Pc], later an independent and more general proof proof was given [ESTV] with no lower 
bound on the window width. The present case is more complicated because the effective coupling 
strength a — ao can be sign-changing. We shall show that it is sufficient if it is negative in the 
mean creating thus a locally average stronger tunelling between the two channels: 



Theorem 3.1 Assume that (i) a — ao £ Lj oc (M), (ii) a(x) — ao — 0(\x\~~ 1 ~ £ ) for some e > as 
\x\ — ► oo. If j R (a(x) — a )dx < 0, then H a has at least one isolated eigenvalue below its essential 
spectrum. 



Proof: We use a variational argument whose idea comes back to [GJ|; sec also [DE, RB|, 



and [ESTV, Sec. Ill]) for a coupled waveguide system. First of all, the assumption (ii) tells us that 



limui^oo \x\ +e (a(x) — ao) = 0, i.e., to any 6 > there is as > 1 such that 



\x\ > as \a(x)—ao\ < 



|l+e ' 



(3.1) 



It is useful to introduce a shifted energy form: for an arbitrary ^ S D(q a ) we put 

Q*[*] :=q a [m-vi(ao)\\nh 



(3.2) 



since the essential spectrum of H a starts at vi(ao), we have to find a trial function "J such that 
Qa^] is negative. We obtain it by a suitable modification of the function (ar, y) := Xi(ui a o) 
which formally annuls (3.2) for a = ao but does not belong to I? . The trial function has to decay; 
in order to make the positive contribution from its tails to the kinetic energy small, we employ an 
exterior scaling. We choose an interval A := [—a, a] for some a > 1 and a function <p £ 5(M) in 
such a way that ip(x) < 1 and f(x) = 1 on A. Then we can define the family {ip a : a £ R} by a 
scaling exterior to A: 



(Pcr(x) 



ip(x) 



if |x| < a 



ip(±a + o~(x =F a)) if ±i > a 



(3.3) 
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By construction, ^^(x)! < 1 holds for all x € R. The sought trial function will be chosen in the 
form ^>(x,y) :— (p (T (x)xi(yj < ^o)- We employ the relations H^vlll = "H^IHj an( i 

= q ao m + / (a(x)-a Q )\^( Xl 0)\ 2 dx, 
Jw 



9« [*] = llvvlli + ^i(ao)||y CT || 
the last one of which is obtained by tedious but straightforward calculation. This yields 

Q«[*] = ^IMll + |xi(0;« )| 2 / (a(:r)-ao)ba^)| 2 ^- (3.4) 



We split now the integration region into two mutually disjoint parts, A and R\.4. Using ( |3.l| ) 
together with the above mentioned bound on ip a we arrive at the estimate 

Q a [9] < a\\<p\\l + 4<5|Xl(0 ; a ° )|2 + |xi(0; «o)| 2 / (a(x)-a ) dx . (3.5) 
ea£ Jr 

By assumption we have J" R (a(x)— ao) dx < and since Xi(0; a o) is nonzero, the last term is negative; 
it is then enough to choose 8 and a sufficiently small to make Q a [\&] negative. | 

Remark 3.2 A case of particular interest concerns weakly coupled Hamiltonian of the type (|j~l|), 
i.e. the situation when a differs from ao only slightly. In that case one can develop a Birman- 
Schwinger analysis in order to derive the perturbative expansion of the ground state energy in 
terms of a parameter measuring the "smallness" of a — a®. This will be done in a separate paper 



[EK]; here we just borrow a result for a further use in this work. 

There are different ways in which a — ao can be small. Suppose that the support of the 
perturbation shrinks, i.e. introduce a a {x) :— a{x/a) with the scaling parameter a 6 (0,1] and 



consider the limit a — > 0+. We have the following result [EK 



Theorem 3.3 Suppose that a— ao is non-zero and belongs to L 1+£ (R, dx) HL 1 (R, |x| 2 c?a;) for some 
s > 0. Then H atT has for small a at most one simple eigenvalue E(a) < v\{ao), and this happens 
if and only if J (a(x) — ao)dx < 0. If this condition holds the following expansion is valid 



y/v x - E(a) = -~ |xi(0; a )| 2 / (a(x)-a )dx 

—r lxi(0;ao)| 2 V IXn(0;a )| 2 / (a(x)-a ) (a(x')-a ) dxdx 

4 fri Jr 2 V v n-vi 



n=2 

T 3 



+0{a i ). (3.6) 



4 A "rectangular well" example 

To illustrate the above result and to analyze the behaviour of coupled waveguides in more details 
we shall now investigate an example. We choose the barrier function a so that the corresponding 
Schrodinger equation can be solved numerically; this happens if a is a step-like function which 
makes it possible to employ the mode-matching method. The simplest nontrivial case concerns a 
"rectangular well" of a width 2a > 0, 



a(x) 



ct\ if |x| < a 
ao if \x\ < a 



for some ai, ao S R. In view of Theorem 3.1 this waveguide system has bound states if and only 



if ai < ao- In particular, one expects that in the case when ai = and ao is large positive the 
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spectral properties will be similar to those of the impenetrable barrier situation studied in [ ESTV | . 
On the other hand, the mode-matching method allows us to treat on same footing the scattering 
processes in our waveguide. Then there is no need to impose the above condition, because the 
"barrier" situation, a\ > ao is expected to exhibit nontrivial scattering behaviour as well. 

Henceforth, we shall denote the transverse eigenvalues in the two regions as v s n := v n (a s ), 
s := 0,1, n € N\{0}. In view of the natural mirror symmetry with respect to the y-axis we may 
consider separately the symmetric and antisymmetric solutions, i.e. to analyze the halfstrip with 
the Neumann or Dirichlct boundary condition at x — 0, respectively. For the sake of simplicity we 
shall also restrict our attention to the case min{ao7 a i} > ®-m '■= + d^ 1 ), when the lowest 

transverse eigenvalue is positive everywhere in the waveguide. The considerations presented below 
remain valid even without this assumption; one has just to replace the trigonometric ground-state 
eigenfunction for hyperbolic which makes the formulae cumbersome. 

4.1 Bound states 

Let us first derive an estimate which allows to localize roughly the eigenvalues. It is based on 
a bracketing argument similar to that used to specify the essential spectrum at the beginning of 
Section p[ The Hamiltonian can be squeezed between a pair of operators, H a N ^ < H a < Ha , with 
additional Dirichlet/Neumann "cuts" at segments perpendicular to the waveguide axis, x = ±o. 
The spectra of the estimating operators can be easily found and sought estimate comes from the 
eigenvalues of the middle part situated below in combination with the minimax principle. In 
particular, we find that the number N of isolated eigenvalues satisfies the bounds 



N D + 1 > N > N D := 



2a 



where [•] denotes the entire part; this complements Theorem 3.1. Furthermore, the n-th eigenvalue 
E n of H a is estimated by 



while the critical halfwidth a n at which the n-th eigenvalue emerges from the continuum satisfies 
the bounds 

(n — 1)tt nn , , . 

. = <a n < 4.2 

After this preliminary, let us pass to the mode-matching method. We start with the simpler case 
when the waveguide exhibits a mirror symmetry w.r.t. the x-axis, i.e. d\ = g?2 = d. 

4.1.1 The symmetric case 

If v = 1, the Hamiltonian decouples into an orthogonal sum of the even and the odd parts, the 



spectrum of the latter being clearly trivial — cf. Remark 2.1. The two symmetries allow us to 
restrict ourselves to the part of SI in the first quadrant, with Neumann or Dirichlet condition in 
the segment (0, d) of the y-axis, and take the transverse eigenvalues determined by the spectral 
conditions 

— ao — 21 cot id if x > a 
-ai = 21 cot Id if < x < a. 

The corresponding transverse eigenfunctions are 

Xn ■= -N" sin y/v°(y-d) if x > a, 

K := ~N^sm^{y~-d) if < x < a, (4.3) 



G 



where N% is a normalization factor chosen to make Xn-, 



(KY 



unit vectors in L 2 (0,d), i.e. 



2^/v^d — sin 

The overlap integrals of elements of the two bases are easily seen to be 



N° N 1 



v\ sin \ d cos \/ ii d 



sin \/ v\ d cos \ v: 



(4.4) 



(4.5) 



A natural Ansatz for the solution of an energy EG v®) is 

*./«(*, V) = E 6H /a e-«-( x - a) Xn(tf) for 



®a/a(x,V) = E a « 



/a 



cosh p^ff 
cosh p w a 

sinh p n a; 
sinh p n a 



(4.6) 



<Pn(y) 



for < x < a 



where the subscripts (we will omit them for the most part) s, a distinguish the symmetric and 
antisymmetric case, respectively. The longitudinal momenta are defined by 



Pn 



As an element of the domain (2.3), the function W should be continuous together with its normal 
derivative at the segment dividing the two regions, x = a. Using the orthonormality of {xn} we 
get from the requirement of continuity 



In the same way, the normal-derivative continuity at x = a yields 



OO S 



tanh 
coth 



Substituting from (4.7) to (Eq), we can rewrite it as an operator equation 

Ca = 0, 

where 



q„ 



Pn 



tanh 
coth 



(PnO) {Xm^n) 



(4.7) 

(4.8) 

(4.9) 
(4.10) 



with the overlap integrals given by (4.5). 

It is straightforward to com pute the norms of the functions (4.6); since n~ 1 q n an n~ 1 p n tend 
to ^ as 7i — ► 00 (see Lemma 2J2 1.), the square integrability of ^ requires the sequences {a n } and 
{b n } to belong to the space £ 2 (n~ 1 ). 



To make sure that the equation (|L9|) makes sense, it is enough to notice that if \& is an 
eigenvector of H a , it must belong to the domain of any integer power of this operator. It is easy 
to check that 

Vi G N\{0} : * G D(H l a ) {a„}, {&„} G £ 2 (n 4i_1 ) ! ( 4 - n ) 

hence the sought sequences should belong to £ 2 (n r ) for all r > —1, i.e. both sequences have a 
faster than powerlike decay. This fact also justifies a posteriori the interchange of summation and 
differentiation we have made in the matching procedure. Furthermore, one can use it to check 
the existence of a convergent series of cut-off approximants to the solutions in the same way as 
in pSTVl, Sec. IV. 1]. 
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Remark 4.1 (an alternative me thod) We can use the orthonormality of {4>n} instead of {xn} 
and express {a n } in analogy to (4/7), and then substitute it into (18). We find that the coefficient 
sequence {b n } is then determined by the following equation: 



b + Kb = . 



(4.12) 



where 



K, 



1 °° f 

n ■= y^(Xm, $k) Pk \ 

qm fc=l L 



tanh 
coth 



(Pfcffl) {4>k,4>n) ■ 



(4.13) 



The two approaches (4.9), ( 4.12| ) are, of course, equivalent, however, it may be useful to combine 
them in order to get a good idea about the numerical stability of the solution. For instance, in 



the situation of [ESTV| the approximants of (4.9) approach the limiting values from above, while 
those referring to (4.12) are increasing. 



4.1.2 The Asymmetric Case 



Let us pass now to the case, when the widths of the ducts are nonequal, v =/= I. In view of the 
mirror symmetry, we shall consider right-halfplane part of Q only with the Neumann and Dirichlet 
condition on the segment [~d,2,di] of the y-axis. The asymmetric case differs from the previous 
one just in the choice of transverse basis: now we can take 



Xn(y) 
<i>n{y) 



Xn(y;a ) 
Xn(y;a 1 ) 



if x > a 

if < x < a, 



where Xn(',(* s ) are of the form (2.7) with the norms 2V* given by (2 
eigenvalues i/„ , v\ are then determined by 



(4.14) 

The corresponding 



-cto = I (cot id\ + cot ld<i) 
-otx = I (cotld\ + cot&fe) 



if x > a 

if < x < a . 



(sec (2.6)) and the overlap integrals are 



N° N 1 / , — , — , — , — 

m n f I ' .1 „■ / O J - /nj • /l 
"0 1 \y v n SinvC»l SHl vC"2 81n V^n 



D 



sin J uld±, sin J vld 2 sin J v®D 



(4.15) 



The res t of the argument does not change and one has again to solve the equation (4J)) (respec- 
tively, ( 4.12 )). By a straightforward modification of the above argument, one can also check that 
the coefficient sequences have a faster than powerlike decay and that the equation can be solved 
by a sequence of truncations. 



4.2 Scattering 

As we said in the opening of this section, the scattering can be treated in an analogous way. 
The incident wave is supposed to be of the form e~ lkrX Xr{y', a-a) , «-e., to come from the right 
in the r-th transverse mode; we have introduced the effective momentum k r := \/k 2 — ^{?. We 
denote by r rn ,t rn , respectively, the corresponding reflection and transmission amplitudes to the 
n-th transverse mode. Due to the mirror symmetry, we can again separate the symmetric and 
antisymmetric situation w.r.t. x — and to write 

r rn = ^(Prn Prn) J ^rn = 7^ (Prn ~ Prn) 1 (4.16) 
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Figure 2: Bound state energies vs. the halfwidth a in the symmetric case for a — 10 5 (□), 50 (o), 
10 (*), 5 (x), 2 (+), 0.5 (.). 



where p° n , u = s,a, are the reflection amplitudes in a half of our waveguide with the Neumann 
and Dirichlet condition at x — 0, respectively. We use the following Ansatz 



*./«,(*, y) ^Ef^-^+^V'-^jxnfe) for x>a 



cos p n x 
n cos p n a 

^ s /a(x,y) =E a » )<t>n{y) 

n=l J sinp^a; 

sin p n a 



(4.17) 



for < x < a 



for the total energy k 2 of the incident wave in the r-th mode. The quantities p n :— \Jk 2 — v\ are 
effective momenta in the "interaction" region. Matching these functions smoothly at x = a we 
arrive in the same way as above at the equation 



Ca = f. 



(4.18) 



where 



C„ 



'/.',» -f Pn { _^ t } (p n a) ) (x m-j fin) 



(4.19) 
(4.20) 




Figure 3: Ground-state eigenfunctions in the symmetric case for a/d — 0.15. 



where the index r corresponds to the incident wave and the overlap integrals are given again 
by (4.1S). The reflection amplitudes are then given by 

oo 

Prm = — $rm + / ] <hi{Xm, <j>n) 5 (4-21) 



they determine the full S-matrix via ( 4.16| ). A quantity of direct physical interest is rather the 
conductivity given by the Landauer formula. If we express it in the standard units 2e 2 /h, it equals 



[k] 



E T^\trnn(k)\ 



m,n—l 



(4.22) 



where t mn (k) are the coefficients (4.16). The summation runs over all open channels. Another 
physically interesting quantity is the probability flow distribution associated with the generalized 
eigenvector = ^ s + ^a, which is defined in the standard way, 



j(x) = 2Im(*(x)V*(f)) 



(4.23) 
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10 



12 



In a„ 



Figure 4: Narrow- window asymptotic power and coefficient as functions of ao- 



4.3 Numerical results 

Since the spectrum behaves naturally at scaling transformations it is reasonable to solve the equa- 



tions ( |4.9| ) and (4.18) in the natural non-dimensional quantities. We mark them by tilde and use 
them to label the axis in the figures, e.g., a = da, a s = a s /d, E = (ir/d) 2 E or k = (ir/d)k. 



4.3.1 Bound states 

Eigenvalues: Figure || shows the bound-state energies as functions of the "window" halfwidth 
a for an "empty window" , oil = 0. Several curves referring to different values of the barrier 
coupling constant ao are plotted. In accordance with the general results of Section 4.1 the energies 
decrease monotonously with the increasing "window" width and one can sandwich them between 
the estimates (4.1). We also see that for a fixed a the energies increase with respect to ao and v; 



recall that their number increases as a function of ao but it decreases as the waveguide becomes 
more asymmetric — these facts are clear from (4.1), (4.2), and LemmaO. It is illustrative 
to confront our results for large ao with the energies computed in [ESTV] for the case which 
corresponds formally to ao = oo. Comparing Figure || with Fig. 2 of the mentioned paper we see 
that our result for ao = 10 5 is practically identical with the latter. 



Eigenfunctions: The evolution of the ground-state wavefunction with respect to ao for an empty 
window of the fixed halfwidth a = 0.15 is illustrated on Figure |3|. If the barrier tunneling is 
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Figure 5: The eigenfunction of the second excited state for v = 1, a/d = 5, a% = 0. 



negligible, a® = 10 5 , the picture is indistinguishable from Fig. 3 in [ESTV]. As ao becomes smaller 
we see how the wavefunction part penetrating the barrier grows. 



Threshold behaviour: Consider again the empty-window case, ot\ = 0. As a consequence 



of Theorem 3.3 we get for the ground-state energy for any fixed ao and a narrow window the 
asymptotic formula (c/. (|3.6|)) 



E(a) 



0(a 3 ) 



c := a 2 al |xi(0;a o ) 



(4.24) 



On the ot her hand, in the case of a window in the Dirichlet barrier, ao — oo, it was conjectured 
in [ESTV] that we may suppose 



E{a) = G) 2 ~ CHfl4 



0(a 5 



(4.25) 



as a — > 0+. The conjecture is supported by a two-sided asymptotic estimate [ EV2] |: there are 
positive ci , C2 such that 



-cia 4 < E(a) 



< 



-c 2 a 



(for a generalization of this resu lt to a larger number of windo ws and higher dimensions see |EV3|). 



Quite recently, a proof of ( 4.25 ) has been proposed by Popov Pop 



This seems to be a paradox. In order to make sense of these considerations, we suppose that 
E(a) = — ca@ for small a of an interval 0.016 < a < a max (a max = a max {a.Q) is chosen in such 
a way to include the best correlated points), and investigate numerically the dependence of the 
coefficients (3 and c (c = d^ +2 c/ir 2 ) on a®. The power like asymptotic behaviour is confirmed when 
we redraw the first eigenvalue curves of Figure ^| in the logarithmic scale. The obtained dependence 
of the coefficients on ao in the symmetric case v = 1, is illustrated on Figure |]. We see that the 
power reaches the values (3 — 2,4 for small and large ao, respectively (a slight shift in the first 
graph is due the truncation; the convergence becomes very slow for small a). At the same time, 
the numerically found c for small ao coincides with that of (4.24). 
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Figure 6: The maximum bending of nodal lines of third eigenfunctions in the symmetric case as a 
function of ao for a fixed window width. 



Of course, the asymptotical behaviour is governed by ( 4.24 ) for any finite ao. The above 
result says only that the transition from biquadratic to quadratic asymptotics occurs for large ao 
at values of a still smaller than those we have used. 



Nodal Lines: In Figure || we plot the third eigenfunction. Its nodal lines are almost straight 
showing thus that the "spikes" at the window edges act almost as a hard barrier. On the other 
hand, a simple argument based on the reflection principle shows that the nodal lines cannot be 
straight. Closer inspection shows that they have the form of a bow bent outward. The maximum 
bending is shown on Figure [| It decreases rapidly with the window width which confirms the 
tunneling nature of the effect. Nodal lines of higher eigenfunctions exhibit (as functions of ao) 
irregularities connected with the changes in the number of the modes. 



4.3.2 Scattering 

Conductivity: Figure ^ illustrates the evolution of the conductivity for the particle coming from 
the right and leaving to the left as a function of the momentum k and the width di- We see that 
the perturbation, — ao in the window, deforms the ideal steplike shape with jumps at transverse 



thresholds; deep resonances are clearly visible, 
practically reproduce Fig. 5a of [ESTV|. 



For an almost impenetrable barrier, ao = 10 , we 
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o =1 00000 




Figure 7: Right-left conductivity as a function of k, di for d — w, a = 2ir, a± = 0. 



Probability flow: Examples of the quantum probability flow are shown on Figure The flow 
patterns change with the momentum of the incident particle and the value of cxq. They exhibit 
conspicuous vortices at the resonance energies which correspond to the "trapped part" of the 
wavefunction. An interesting phenomenon is illustrated on the first two graphs of Figure [| for 
cto — 10 5 there is a double vortex (corresponding to the sharp stopping resonance of Figure |?|), 
right-handed in the upper duct, while for ao = 50 we get a left-handed vortex. The conductivity 
is small in these situations so the waveguide system is effectively closed for the particle transport. 
As ao decreases the conductivity grows and the waveguide opens — cf. Figure H for ao = 10, 2. 



14 



oe =1 00000 



t t ( s -* — 



Iff / , — -~*\v » i ' 




2 
1 

y o 
-1 

-2 



a Q =50 




i iff > — _ — *^ v * 

; i J j J / f r ~ > > > » 




a =10 





y o 



-3 



-2 



a Q =2 



-2 



Figure 8: Probability flow patterns for fc = 1.745 and ai = in the symmetric situation. 



5 An Aharonov-Bohm cylinder 

In the closing section we want to show now how the preceding considerations modify for a different 
geometry: we consider a nonrelativistic quantum particle living on the surface of an infinite straight 
cylinder of a radius R, which is threaded by a homogeneous magnetic field B parallel to the cylinder 
axis. We assume that the motion is further restricted by a (5-barrier supported by a line parallel 
to the axis. 

5.1 General considerations 

The configuration space is sketched on Figure ^ where we indicate also how the coordinate system 
is chosen. In these coordinates we have 

n:={(x,y,z)£R 3 \ y 2 + z 2 = R 2 } . (5.1) 

Choosing the gauge so that the electromagnetic potentials fulfill <p(x) = and A(x) = \B x x, the 
Hamiltonian acquires the form 

H a := HV + A) 2 = -A - i^(yd z - zd y ) + (5.2) 
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B 




Figure 9: Cylindrical strip with a 8 barrier in axial magnetic field. 



away of the barrier, where B := \B\; as before we put H = 2m = 1, and also e = —1 having in mind 
an electron. The subscript a indicates the real function which defines the shape of the barrier 



as in the strip waveguide situation; it will enter the boundary condition (5.4) below. The vector 
potential has obviously the the angular component A v = A and it equals 



(5.3) 



where 4> is the magnetic flux through the cylinder. Recall that in the rational units we use here, the 
natural unit of the magnetic flux is <po := (2tt)^ 1 . Since we deal with a quantum system living on a 
surface it is natural to "unfold" it and to study (5.2) on a planar strip with appropriate boundary 
conditions, namely 



ip(u,0+) = ip(u,2irR-) =: ip(u,0), ip v (u,0+) - ip v (u,2wR-) = a(u)ijj(u,0), 



(5.4) 



where the subscript denotes again a partial derivative. To this aim, we introduce the unitary 
transformation U : L 2 (&) — > L 2 (R x [0, 2nR],dudv) by 



v v 
(Uip)(u,v) := ip(u, Rcos — , i?sin — ) , 
R R 



(5.5) 



which maps f2 onto the strip f2 



: [0, 2itR\; the operator H a is then unitarily equivalent to 



H a := UH a U- 1 = -d 2 u + {-id v + A) 2 , 

with the domain 

D{H a ) := {i> G W|(n) | Vu G E : b.c. (5.4) are satisfied } . 
We will need also the quadratic form q a associated with H a . Its domain is D(q a ) 



Q 

2iA 



\Vip\ 2 (u,v)dudv + / a(u)\ip(u,0)\ 2 du 



(ipd v ip)(u,v)dudv + A / (u,v)dudv. 
n Jn 



(5.6) 

(5.7) 
Wi(Q) and 

(5.8) 



As a comparison operator we employ again the one with a(u) = a = const when we can solve 
the Schrodinger equation by separation of variables. We denote by {t^l^Li and {x n }^Li the 
(properly ordered) sequences of the transverse eigenvalues and the corresponding eigenfunctions, 
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respectively. Since our system is now more complicated due to the presence of the magnetic field, 
we have to distinguish several possibilities: 



1. No barrier, a = 0: 



Wa: xe(v) 



1 



(5.9) 



The tilde marks the eigenfunctions corresponding to the eigenvalues (-^ + A) 2 ; to get {v n } one 
has to arrange thee latter into a ascending sequence. The respective eigenfunctions will be then 
denoted as {Xn}- 



2. a/0 and 2RA g N: 

VnGN\{0}: X n(v) = N n e 
where N n denotes the normalization factor chosen to make \n a unit vector in L 2 (0, 2irR), 



\N n \ 2 := 



~ [l-cos27rJ?.( v ^vT+ J 4)] 



47r_R v /F^"(l-cos 2ir_RAcos27ri?, v /i^r)+2sin 27r.R v /7^"(cos 2-7T.R.4-COS 

and the increasing sequence {v n } arises from the spectral condition 

cos 27r i?yT 



a = 2l\ cot 2nRl 



sin 2nR£ 



(5.10) 



(5.11) 



(5.12) 



In analogy with Lemma |2~^ (a) we have ^/u^ £ ^n( n ~ 1; n ) f° r an y 71 G N\{0}. 
3. a ^ and integer flux , 2RA £ 2 N: 

The transverse eigenfunctions are of the form (5.10), while the spectral condition ( 5.1 2| ) changes 
to 

a = 2tt&i\%Rl. (5.13) 
Moreover, for I £ 4 N we always get the trivial solutions 



xT v (v) 



-iAv 



'■kR 



sin Iv , 



(5.14) 



which are independent of a. The roots {vn} of (5.13) satisfy the estimates \fv\ £ ^(0,1) and 
e 27?(2n-3,2n-l) for n £ N\{0, 1}. 

4. a/0 and half-integer flux , 2RA £ 2N+1: 

As in the two preceding cases the transverse eigenfunctions are still (|5.10 ) but the spectral condi- 
tion ( 5.12 ) changes now to 

- a = 2£cot7ri?£. (5.15) 



The wave functions of the trivial solutions, t £ t^(2N+1) are (5-14) again and the nontrivial roots 
of ( |5.15| ) satisfy y/^£ ^(2n-2,2n) for all n £ N\{0}. 



Remark 5.1 The integer and half-integer values here refer to the natural flux unit mentioned 
above. The essential instrument for proving the existence of bound states is the requirement 
Xi(0) 7^ (compare, e.g., to Eq. ( |3.4[ ) above). In the absence of a barrier, a = 0, the eigenfunctions 
are always positive (see (|5.9| )) and Xi(Q) ^ holds for a < a m := — 2 sm v ^ RA ■ In the case of a 
(half-)integer flux we have to ex clud e the trivial solutions ( 5.14 ). It is an analogy of the trivial-part 
exclusion described in Remark |2.l| ; the difference is that the triviality now does not come from 
the waveguide geometry, but rather from the magnetic field, i.e., an external parameter. It is also 
clear that the ground-state eigenfunction of the class ( 5.10] ) can vanish at the barrier only if a > 
and the flux is half-integer. 
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Figure 10: The quantum probability flow on the cylinder surface for A = 0.5, k = 1.705 and 
fii = 10" 5 . 
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5.2 An example 



We shall illustrate the above consideration on a "rectangular-well" example analogous to that of 
Section |J. Most of the argument proceeds as there, one has just to use the different eigenf unctions 
and to recompute for them the overlap integrals: 



sin 2irRy/v^ (cos 2irRA — cos 2-KR\pt^ l ) 
v\ sin 2-kR^/v^ (cos 2irRA - cos 2nRy^ ) 



for a ^ ^ ai with m, n £ N \{0}, and 



(Xm,<j>t) 



4N° 



sin 2-kRA + i cos 2irR 



(5.16) 



(5.17) 



for a 7^ 0, a% = with m S N\{0}, <£Z- Note that in the latter case one has to substitute the 
ordered basis {4>n}%Li (together with the corresponding eigenvalues) into Ansatz ( jOj ) to make the 
numerical procedure of cut-off approximations convergent. 

Unless «i < we have to exclude here the possibility 2RA £ 2 psf+1 again (c/., e.g., (fh^)). Next 
we can restrict our atte ntion on ly to the situation when 2RA E [0, 1 ) U ( 1, 2) because A appears 
in the overlap integrals ( 5.16|) , ( 5.17 ) and in the spectral condition ( 5.12 ) as an argument of the 
periodic functions sin, cos; the integ rals and the transverse eigenvalues are the only quantities 
which affect the equations (19) and ( 4.18 ). 

I n fac t, we can take 2RA from [ 0,1) only because the replacement 2RA i— > 2 — 2RA in ( 5.16| ) 
and ( 5.17 ) (the spectral condition ( [5.12 ) does not change at all) is equivalent to the exchange 
A i ► — A which coincides with the conjugation of Hamiltonian (|ij). It is well known fact that 
such an operator has the same energies while the corresponding eigenfunctions are given by a 
simple conjugation. 

As an example, the evolution of the quantum probability flow w.r.t. c<q is illustrated on 
Figure |Tc|. 
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Ground-state cigenfunctions in the symmetric case for a/d — 0.15. 
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Probability flow patterns for k = 1.745 and ct\ = in the symmetric situation. 
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Cylindrical strip with a S barrier in axial magnetic fieldj 
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The quantum probability flow on the cylinder surface for A — 0.5, k = 1.705 and 
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